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We apply and test the recently proposed "extended scaling" scheme in an analy- 



Abstract 
> 

■^j- | sis of the magnetic susceptibility of Ising systems above the upper critical dimension. 

The data are obtained by Monte Carlo simulations using both the conventional Wolff 
cluster algorithm and the Prokof 'ev-Svistunov worm algorithm. As already observed 
for other models, extended scaling is shown to extend the high-temperature critical 
qq ' scaling regime over a range of temperatures much wider than that achieved conven- 

^^ . tionally. It allows for an accurate determination of leading and sub-leading scaling 

indices, critical temperatures and amplitudes of the confluent corrections. 
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1 Introduction 

An approach to critical-point scaling motivated by high-temperature series expansions 
has recently been developed, which aims to extend the scaling window well beyond the 
critical regime pQ. At a temperature T sufficiently close to the critical temperature T c , 
divergent thermodynamical averages display the scaling behavior 0(T) ps Aot~ p where 
Ao is a constant amplitude, p is the critical index associated with the observable O and 
t = (T — T c )/T c is the standard reduced temperature. In order to improve the temperature 
range where scaling holds, other thermal scaling variables have been considered in the 
literature. In particular the alternative reduced temperature 

T= T-T 1 = 1 -l 
T & 

is popular in analysis of series expansions. This variable is also natural in renormalization 
group analyses where the temperature variable usually appears through j3J (J being a 
coupling constant between microscopic degrees of freedom), and thus the reduced variable 
is defined as in (11.11) . The relationships between the two reduced temperatures are 

t = -^ = r + r 2 + r 3 + ... , (1.2) 

1 — T 

r = J— = t-t 2 + t 3 -..., (1.3) 

such that in the vicinity of the critical point 

0{T)^A t- p or 0{T)^A T- p . (1.4) 

While both t and r vanish as the critical point is approached, they have very different 
high-temperature limits, 

lim t — oo , (1.5) 

T — >oo 

lim r = 1 . (1-6) 

T — >oo 

In [T], a rationale was given explaining why the alternative reduced temperature r in 
(11.11) may be superior to the more traditional variable t in scaling analyses, at least in the 
high-temperature regime. This rationale stems from the following observations. Firstly, 
suppose the constant A Q in (11.41) is promoted to a temperature-dependent amplitude as 
follows, 

0{T) oc 0*{T) ~ T^° (T - T c )~ p ^ T^°- p (l -^] " ~ ^°t~ p , (1.7) 



where 0o = P — 4*0 ■ From ( 11.61) . 



T 



lim 0*{T) ~ (3^° , (U 

T—>oo 



and 0o may be chosen so that 0*(T) matches the high-temperature series expansion 
(HTSE) for 0(T) in this limit. Note that a similar approach cannot be implemented for 
scaling in t. I.e., scaling in r, unlike scaling in t, allows 0*(T) to represent the correct 
asymptotic behavior of 0(T) in the T — >• oo limit as well as close to criticality. Inspired 
by fll.7p we now write the full expression for the observable 0(T) as 

0{T) = A 0*(T) = A ^ T~ P + ■ ■ ■ , (1.9) 

where the dots represent higher-order additive corrections. 

In [H |5], Campbell et al. proposed and tested the extended scaling scenario in the 
Ising model in two dimensions as well as the Ising, Heisenberg and XY models in d = 3 
dimensions. These works convincingly established the superiority of the method over the 
conventional scheme between the lower and upper critical dimensions. The inherent ideal- 
ness of the extended scaling approach to study systems below the lower critical dimension 
was demonstrated in [3J. In this work, we test extended scaling for the <i-dimensional 
Ising models with d = 5, 6, 7 and 8, i.e. above the upper critical dimension d c = 4. Fur- 
thermore, while the d — > oo (mean-field) case was also investigated in [2] (where extended 
scaling is exact for all T > T c ), here we include a detailed analysis of corrections in the 
extended scaling scheme. 

In the next section we review the known results about the leading and sub-leading 
scaling behavior of the Ising model above d c . We then define the extended scaling of 
this model, taking into account scaling corrections. In Section 3 the MC approach is 
presented and in Section 4 the resulting data is discussed. Agreement with results from 
HTSE confirms the efficacy of extended scaling, which is then used to determine the 
critical parameters governing scaling in higher dimensions. Since above d = 6 dimensions, 
the dominant corrections to scaling are due to analytic terms, this work provides a full 
account of leading and confluent corrections in the susceptibility above the upper critical 
dimension. We conclude in Section 5. 



2 Extended scaling in high dimensions 

Above the upper critical dimension d c = 4, the leading scaling behavior for the magnetic 
susceptibility is given by mean field theory and is characterized by the exponent 7 = 1. 
The additive correction-to-scaling terms are expected to be non-classical. The generic 
scaling form for the reduced susceptibility for the Ising model is 

x(T)/(3 = rr' y (1 + M e + b 2 t 2e + ... + Cl t + c 2 t 2 + ...), (2.1) 

where the dots indicate higher-order terms. In (12.11) . the confluent (or non-analytic) 
corrections involve the universal exponent 6, while the remaining correction terms are 
analytic. (While in principle there could be several confluent correction exponents in- 
volved in (12.11) . it is sufficient in this work to consider the single exponent 9.) Analysis 
of the </> 4 model above d c = 4 shows that the dominant critical behavior is determined by 
the Gaussian fixed point, leading to \ ~ (0 2 ) ~ ^ _1 an d £ ~ t' 1 ' 2 . The confluent scaling 



Table 1: Selection of previous estimates for the critical temperature T c for high- 
dimensional Ising models together with the refined extended scaling estimates obtained 
in Sec. 4 of this paper. Here, FSS means a Monte Carlo approach using finite-size scaling. 



Reference method & year 


T c 


T c 


T c 


T c 




d = 5 


d = 6 


d = 7 


d = 8 


[H] l/d expansion (1964) 


8.7881... 


10.8397. . . 


12.8712... 


14.8923. . . 


[I] (critical expansion) (1981) 


8.8162(5) 


10.8656(8) 






[7] HTSE (1993) 


8.777(2) 








[6] HTSE (1993) 




10.8348(4) 


12.8690(4) 




[9] FSS (1999) 


8.77844(2) 








[10] FSS (2000) 




10.835(5) 






HI] FSS (2001) 






12.870(5) 




[12] FSS (2005) 


8.7785(5) 








[13] FSS (2006) 






12.870(5) 


14.893(3) 


This paper (2008) 


8.7777(9) 


10.8353(4) 


12.8690(3) 


14.8933(8) 



4 term. The perturbation felt by the system inside 
l ~ 2 ) 2 V£ (since the average is taken over a Gaussian 

.-2+d/2 



corrections are due to the irrelevant 

the correlation volume is (0 4 )/£ ~ 

distribution), which thus behaves in the vicinity of the critical point as x 2 /£, d 

As a consequence, the exponent 9 is given by 



e 



d-A 



(2.2) 



for d > d c — 4. Thus, for d > 6, where 9 > 1, the confluent corrections are expected to 
be overwhelmed by the analytic ones. In six dimensions, the corrections are expected to 
be modified by a logarithmic factor [4J [5] , and, explicitly, 



X(T)/P = 


Tr l + Bt~2 +C + Dt* + ... , ford = 5, 


(2.3) 


x{T)/(3 = 


Tr 1 + B\nt + C + Dt\nt + ... , ford = 6, 


(2.4) 


x(T)/(3 = 


rr 1 + C + Dt* + ..., for d = 7, 


(2.5) 


x{T)/(3 = 


rt" 1 + C + Dt + . . . , for d = 8. 


(2.6) 



Series analysis techniques have been used to verify the leading correction-to-scaling 
exponent 9 in five and six dimensions and to determine some of the critical amplitudes 
[4]. In particular, accepting the classical value 7 = 1, Guttmann measured 9 = 0.50(5), 
T = 1.311(9) and B = -0.48(3) in five dimensions as well as 9 = 0.98(5) and T = 1.168(8) 
in six dimensions [1]. It was not possible to determine the subdominant amplitudes, nor 
to confirm the logarithmic term, in six dimensions using these techniques [4] . 

The high-temperature series expansion (HTSE) for d-dimensional Ising models was 
given to fifteenth order in tanh(/3) by Gofman et al. in [6], and compared with results 
of Monte Carlo simulations in six and seven dimensions. The HTSE in five dimensions 



was analyzed in a related study [7J. In contrast to five dimensions where the leading 
corrections are confluent, analytic corrections to scaling were favored in d = 6 and d = 7 
dimensions [6J using the HTSE approach. However, the logarithmic term in six dimen- 
sions was also not explicitly handled in the series expansion approach of [Bj. With no 
prior assumptions regarding exponent values, and ignoring the logarithm in six dimen- 
sions, the HTSE gave 9 = 0.45(10), 9 = 1.0(3) and 9 = 0.8(2) in d = 5, d = 6 and d = 7 
dimensions, respectively [6j [TJ. For convenience, a list of previous measurements of the 
critical temperatures in high-dimensional Ising models is given in Table [lj 

Here, we report on simulations of the Ising model in d — 5, 6, 7 and 8 dimensions 
using the new extended scaling method recently developed by Campbell et al. [T]. It is 
claimed that this new method can extend the scaling regime well beyond that traditionally 
available [jQ. This claim has been verified in all models so far tested below their upper 
critical dimensions [2] and the method is particularly well suited to models which order 
at zero temoerature, i.e., below the lower critical dimension [3J. A central theme of this 
paper is to test the applicability of the method for high dimensions, i.e., above the upper 
critical dimension, and if it is proven more suitable than "standard temperature scaling" , 
the method will be used to investigate corrections to scaling and amplitudes beyond the 
leading critical behaviour. 

We consider the magnetic susceptibility x{T)//3 which is defined in the paramagnetic 
phase as 

«=kN. ( 2 - 7 ) 

and numerically estimated by integration of the spin-spin correlation functions over the 
lattice. With complete randomization in the high-temperature limit (T — > oo), (siSj) = Sij 
there, so that x(T)//3 approaches unity. Therefore the HTSE for the reduced susceptibility 
has the form 

X(T)/P = a + oi/9 + a 2 (3 2 + a 3 /3 3 + . . . , (2.8) 

where do = 1 for all d. Comparing this to (11.81) . it is clear that <p x = for the reduced 
susceptibility. This is the reason why the second equation in (11.41) may represent the sus- 
ceptibility over an extended high-temperature range. The limit (II. 5j) , on the other hand, 
restricts the suitability of the conventional variable t for temperature scaling analysis in 
the vicinity of the critical point. Writing the five-dimensional critical expansion (I2.3P in 
terms of the extended scaling reduced temperature r instead of t, one has 

X (T)/(3 = Tr- 1 + Br- 1 2+C(r), (2.9) 

where the higher corrections are contained in the function C(r), which goes to a constant 
C in the high-temperature limit. Since both r and x(T) approach unity there, one has 
C = 1 — T — B. Similar considerations in six and higher dimensions lead to 

X {T)I(3 ~ Ft' 1 + B In (r) + C, (2.10) 

and 

X (T)/i3^Tr- 1 + C, (2.11) 

respectively, where C = 1 — T. 



3 The Monte Carlo algorithms 

The (main) data presented in this work have been obtained using the so-called worm 
algorithm introduced by Prokof ' ev and Svistunov and which provides an efficient Monte 
Carlo sampling of the dominant terms of the HTSE of spin-spin correlation functions [T4"l 
[T5] . The data resulting from this approach are labelled MCHTSE (for Monte Carlo high- 
temperature series expansion) throughout this paper. In the case of the Ising model, this 
HTSE is 

(Sis,) = ^E W^' ^ = 4( c °sh P) dN J2 s i s i II (1 + s^tanh/?). (3.1) 

Here, s; is the spin at site % of the d- dimensional lattice, Z is the partition function, 
the summation is taken over configurations and (k, I) denotes a pair of neighboring sites 
on the lattice. Since J2 Si =±i s i 2p = 2 while Z) Si =±i Si 2p+1 = for integer p, the only- 
graphs that contribute to the expansion are paths joining sites i and j (which we call 
sources) and closed loops. These graphs can be sampled stochastically using the following 
rules: one of the sources, say that initially at site i, is moved to one of the neighboring 
sites, say i', and a bond is added between % and i' if none was hitherto present, or the 
bond is deleted if already present. If the two sources are on the same site, i.e. i — j, 
they are both moved to a new randomly chosen site. These three moves are chosen 
using a Metropolis prescription: the probability to add a bond is min(l, tanh/3J), that to 
delete a bond is min(l, l/tanh/3J) while the probability to move both sources can be set 
freely. This algorithm is known to be slightly more efficient than Swendsen-Wang cluster 
algorithm [16]. To accelerate the dynamics, the state of the bonds starting from the same 
site is stored as a single bit. Furthermore we implemented a continuous-time version of 
the worm algorithm. At each Monte Carlo step, the probability u a of all possible moves 
is calculated. The time r that the system will stay in the same state is evaluated as [IT] 

r=l + Intf— — ^ -) , (3.2) 

where rj is a random variable uniformly distributed over ]0; 1]. The time is increased by 
r and one of the moves a is chosen with probability u a and applied to the system. 

As the temperature approaches the critical temperature, the worm algorithm generates 
graphs contributing to higher and higher orders of the HTSE. The number of these graphs 
grows exponentially fast which may cause critical slowing-down of the convergence of the 
averages. To check the results, additional data have been computed using the standard 
Wolff algorithm [18] whose critical slowing-down is much better understood. For a given 
lattice size, the simulations have been performed at a temperature close to T c , and data 
for nearby temperatures were obtained by the so-called histogram reweighting method 
|19j . The reliability of the MCHTSE data is demonstrated in Fig. (TJ where estimates 
for the susceptibility from this method are compared with those from the HTSE of [6] 
and conventional Monte Carlo data for lattice sizes L = 6-14 in five dimensions. Clearly, 
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Figure 1: MCHTSE estimates for the susceptibility of the five-dimensional Ising model 
compared with conventional Monte Carlo estimates (dot-dashed curves, colored online) 
using the Wolff algorithm together with histogram reweighting and with the HTSE of [6] 
(dashed line, black online). 

the MCHTSE data accurately follow the HTSE curve for large T and are comparable 
to conventional MC estimates closer to T c , where finite-size effects become important. 
Furthermore, the MCHTSE for L > 10 is independent of L down to T rs T c + 0.1 ~ 8.9. 

4 Extended scaling analysis of MC data 

Using the numerical approach presented in Sec.[3l we computed the magnetic susceptibility 
X in Ising models above d c . Extended scaling is tested in Figs. [TH7] and then employed 
to determine the leading corrections for the scaling behavior of the susceptibility. In 
particular, we firstly determine 9, T, B and C in five dimensions. The agreement between 
our measurements there and those of [1] establishes confidence in the extended scaling 
approach. The fifteenth-order HTSE for <i-dimensional Ising models of Gofman et al. in 
[6] is indeed expected to be accurate at sufficiently high temperature, while the critical 
expansion of Guttmann [1] should be reliable close to criticality. Extended Scaling is then 
used to measure 9, T, B and C in six dimensions as well as 9, F and C in seven and eight 
dimensions. The scaling forms fl2.3|) . (j2.4| . (12.51) and (I2.6P are confirmed (including for 
the first time the logarithmic correction in the leading correction-to-scaling term in six 
dimensions) and hence the general formula (j2.2p is supported. These forms are then used 
to determine refined estimates for T c in d = 5 to d = 8 dimensions. The amplitude of the 
logarithm in six dimensions turns out to be small, explaining the difficulties in verifying 
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Figure 2: MCHTSE estimates for the susceptibility of the five-dimensional Ising model 
compared with conventional Monte Carlo estimates (dot-dashed curve, blue online), the 
HTSE of [6] (dashed line, black) and the critical expansion of |3] (dotted line, green online). 
The inset illustrates how the critical expansion becomes inaccurate at high temperature 
where the MCHTSE data cross over to the HTSE curve. 

it numerically almost thirty years ago [3]. 



4.1 Five-dimensional case 

In Fig. [21 the d = 5, L = 14 MCHTSE data are compared with the HTSE 0, the critical 
expansion [I] and conventional L = 14 MC data on a double-logarithmic scale using the 
standard reduced temperature t. The MCHTSE data follow the critical curve for t > 0.05 
to t ^ 0.2, where they cross over to the HTSE curve. The insert in Fig. [2] is a blow-up of 
this region, and clearly illustrates the deviation of the critical expansion from the HTSE 
and the MCHTSE data. 

To extract the parameters governing the scaling in five dimensions, the following ex- 
tended scaling analysis was applied to the MCHTSE data. From (12.31) . the first three 
terms characterizing the critical behavior of the susceptibility are expected to be of the 
form 

X (T)/(3~rr^ + BT» + C, (4.1) 

where $ = —7 + 9 from (12.11) , r is given in (11. ip and extended scaling gives C = 1 — T — B. 
Firstly, for L — 14, a six-parameter fit of the MCHTSE data to (14. ip between T = 8.9 
and T = 25 gives an estimate for T c of 8.7743(95), which compares well with the best 
estimate in the literature (T c = 8.77844(2) [9]). Refined extended scaling estimates for 
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Figure 3: Extended scaling for the susceptibility of the d — 5 Ising model. The extended 
scaling curve (solid, red online) coincides with the critical curve (dotted, green online) [4] 
close to r = (T = T c ) and with the HTSE (dashed, black) [6] close to r = 1 (T — > oo) 
as well as with the MCHTSE estimates in between. 

the critical temperatures at each value of d are given later in this section and listed in 
Tables [1] and El Accepting this more accurate estimate for d — 5, a fit to the remaining 
five parameters yields an estimate for 7 of 1.0(2). Accepting the mean field value 7 = 1 
and fitting to the remaining four parameters gives T = 1.28(2), B = —0.44(19), so that 
1-T-B = 0.16(21), which is compatible with C = 0.17(21). Accepting that C = l-T-B 
(i.e., using extended scaling), a three-parameter fit to the amplitudes T and B, as well 
as to the correction exponent gives d = —0.47(8). Finally, accepting that $ is, in fact 
—0.5, as in (12.31) . a two-parameter fit to the amplitudes of the leading and first correction 
terms gives T = 1.291(3) and B = —0.310(8). The goodness of fit for each of these 
measurements (and for each fit reported below) is monitored through the chi-squared per 
degree of freedom and is observed to be reasonable in each case. These results, together 
with the higher dimensional ones, are summarized in Table [2J where the refined (see 
below) extended scaling estimates for the critical temperature are also given. 

The same analysis applied to the L = 16 MCHTSE data yields very similar results 
and, in particular, T = 1.291(3), B = —0.311(8). For comparison, Guttmann's series 
analysis yielded T = 1.311(9) and B = -0.480(30) [lj. 

This extended scaling (i.e., (EH]) with C{r) = C = 1 -T - B andT = 1.291(3) and 
B = —0.311(8)) is depicted in Fig. [3] for the d — 5 case. This line coincides with the 
HTSE [6] for high temperature (r — > 1) and with the critical expansion [1] close to T c 
(r = 0). Between these extremes, it successfully follows the MCHTSE data. It is clear 
from Fig. [3] that the critical reduced-temperature line of [1] deviates significantly from the 
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Figure 4: MCHTSE estimates for the susceptibility of the six- dimensional Ising model 
compared with estimates coming from conventional Monte Carlo (dot-dashed, blue on- 
line), the HTSE of [6] (dashed, black) and the critical expansion of [lj (dotted, green 
online). The inset illustrates the deviation of the critical expansion from the HTSE at 
high temperature. 

data away from T c , a circumstance typically ascribed to correction terms. Similarly the 
HTSE curve deviates from the data as the temperature is reduced. On the other hand, 
the reduced temperature curve nicely follows the critical expansion for small r and the 
data to very large T. Thus the extended scaling method proposed in [1] is indeed seen to 
be superior to the conventional approach, also in high dimensions. 

In summary, the extended scaling analysis of the MCHTSE data yields similar but 
improved results to that of [I] for the d = 5 Ising susceptibility. Therefore, confidence 
has been established in the method, which can now be applied to higher dimensions to 
determine critical parameters which were unobtainable previously. 



4.2 Six-dimensional case 

In Fig. HI the six-dimensional L = 14 MCHTSE data are compared with the HTSE and 
the critical expansion as well as conventional (reweighted) L = 14 MC data. Similar to 
the five-dimensional case, the MCHTSE data follow the critical curve close to t — and 
switches to the HTSE curve for larger t. The insert of Fig. H] illustrates the deviation of 
the critical curve from the HTSE and the MCHTSE data. 

Although the leading correction term in six dimensions is expected to involve a loga- 
rithm after (|2.4j) . we initially fit to the form (14. ip so as not to a priori bias in favor of the 
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Figure 5: Extended scaling for the susceptibility of the Ising model in d = 6 dimensions. 
The extended scaling curve (solid, red online) traces the critical curve of [4] (dotted, green 
online) close to r = and the HTSE [6] (dashed) close to r = 1 and as with the MCHTSE 
estimates in between. 



14, and fitting to 55 data points 
1]) gives T c = 10.8318(24), 



logarithmic structure. Using the MCHTSE data for L 
from T = 10.88 to T = 30, a six-parameter fit to the form 
close to the value T c = 10.8348(4) measured in [6]. Following the same procedure as in 
the five-dimensional case, and accepting this value, subsequent fits yield 7 = 0.95(25) 
and $ = —0.01(40). Since the latter value is close to zero (and because the fitting pro- 
cess becomes more unstable as the number of parameters is reduced), we conclude that a 
logarithmic term may indeed be present at the leading-correction level. 

We therefore perform a four-parameter fit to the form (12.101) which gives 7 = 0.994(3). 
Accepting the mean field value 7 = 1, a three-parameter fit gives 1 — T = —0.1607(24) and 
C = -0.1605(18). Accepting that C = 1-T (extended scaling) then yields T = 1.1606(17) 
and B = 0.0571(27). The value for T compares well with Guttmann's estimate 1.168(8) 
[3]. The logarithmic term was not addressed in previous analysis [HE]. 

Extended scaling (i.e., fl2TT0|) with T = 1.1606(17), B = 0.0571(27) and C = 1 - T) 
is depicted in Fig. [5] for the d = 6 case. This line coincides with the HTSE [6] for high 
temperature (r — > 1) and with the critical expansion [1] close to T c (r = 0) and traces 
the MCHTSE data in between. 

4.3 Seven-dimensional case 

For d = 7 dimensions, we simulated at three different lattice sizes (L = 10, 12 and 14) to 
again determine the temperature range over which the MCHTSE data are L independent 
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Figure 6: Extended scaling (solid curve, red online) for the susceptibility of the Ising 
model in d — 7 dimensions with the HTSE of [6] (dashed curve, black online) and the 
MCHTSE data. 



(T = 12.89-30). Fitting the L = 12 data to all four parameters in (I2TTT1) yields T c = 
12.8665(5), close to the value 12.8690(4) (obtained in j6] using the HTSE), which we now 
accept (pending our refined extended scaling estimate, which we give below). A three- 
parameter fit then yields 7 = 1.000(2), and accepting the mean field value 7 = 1, we find 
T = 1.1086(9) and C = -0.1215(20). Accepting that C=l-T (extended scaling) finally 
gives T = 1.1008(5). The extended scaling plot for seven dimensions is given in Fig. [61 



4.4 Eight-dimensional case 

A similar analysis in eight dimensions (with L — 8) gives the successive estimates T c = 
14.8893(20) (to be compared with the value T c = 14.893(3) from finite-size scaling [T3]). 
7 = 0.998(2), T = 1.0935(45) and C = -0.0950(26). Extended scaling (C = 1 - T) then 
yields T = 1.0836(5). The extended scaling plot for d — 8 is depicted in Fig. [TJ, where 
it too is compared with the HTSE [6] and the critical expansion [I]. These results, with 
refined estimates for the critical temperature, are summarized in Table |2j 

Having established that 7 = 1 in each case, and the expected scaling forms (12. 3p to 
(12.61) indeed hold, we return to the estimates of T c in each dimension. Refitting to these 
forms for the critical temperature (as well as the amplitudes) yields the more refined 
estimates T c = 8.7777(9), T c = 10.8353(4), T c = 12.8690(3), T c = 14.8933(8), for d = 5, 
d = 6, d = 7 and d = 8 respectively. These values should be compared with previous 
estimates in the literature, which are listed in Table [TJ 

With the Boltzmann constant and the coupling strength having both been set to 
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Figure 7: Extended scaling curve (solid, red online) for the Ising susceptibility in d = 
dimensions with the HTSE of [6] (dashed curve, black online) and the MCHTSE data. 



unity, for a regular lattice of the type considered here the mean field expressions for the 
Ising critical temperature and amplitude are T c = 2d and T = 1. The results for T c 
and T, summarized in Table [2], illustrate this approach to these classical values as the 
dimensionality is increased. The approach to mean field is also illustrated in Fig. [BJ using 
our estimates for T c for d = 5 to d = 8 together with the exact result for d = 2 [20] and 
best estimates from the literature for d = 3 and d = 4 (2TJ [22]. These data are fitted 
to polynomials of increasing degree in 1/d until the goodness of fit becomes reasonable. 
The solid curve represents a fit to the polynomial form /3 C = b\/d-\- b%ld 2 + b^/d 3 + b^/d , 
with 6i = 0.5044(3), b 2 = 0.1935(30), b 3 = 0.3540(101) and b 4 = 1.5335(104). The dashed 
line represents the mean field approximation T c = 2d. 



5 Conclusions 

We have presented estimates for the magnetic susceptibility in the paramagnetic phase of 
high- dimensional Ising models, obtained using the Prokof 'ev-Svistunov worm algorithm 
which enables the stochastic generation of the HTSE by Monte Carlo methods. The 
resulting data are checked against the critical expansion [1], the HTSE [B] and the results 
of a standard Monte Carlo approach using the Wolff algorithm (THJ. When applied to 
this data, the new extended scaling approach of [1] is demonstrated to work well in d = 5 
dimensions, where known results for the estimates of critical parameters [1] are recovered 
and improved. Having thus established confidence in the extended scaling approach above 
the upper critical dimension, it is then applied to the data obtained in d = 6, d = 7 and 
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Table 2: Summary of the numerical results obtained in this paper for dimensions d = 5 
to d = 8. The leading exponent is compatible with the mean field result j — 1. In d = 6 
dimensions, the leading correction is logarithmic. Analytic corrections dominate confluent 
exponents beyond d — 6. 



d 


T c 


7 


•d 


Amplitudes 


5 


8.7777(9) 


1.0(2) 


-0.47(8) 


T = 1.291(3), B = - 
C = 0.019(11) 


-0.310(8) 


6 


10.8353(4) 


0.994(3) 


[logarithm] 


T = 1.1606(17), B-- 
C = -0.1606(17) 


= 0.0571(27) 


7 


12.8690(3) 


1.000(2) 




T = 1.1008(5), C = 


-0.1008(5) 


8 


14.8933(8) 


0.998(2) 




T = 1.0836(5), C = 


-0.0836(5) 



d = 8 dimensions to deliver estimates for critical parameters (especially those governing 
confluent corrections) there. In particular, the logarithmic correction in six dimensions 
(unobtainable in previous analyses [U Ej) is clearly verified. It is also observed how the 
infinite dimensional (d — ► oo) limit leads to the mean field theory. In this way, the 
general formula (j2.2p is supported for the Ising model and a full account of the leading 
and confluent corrections to scaling in the odd sector above the upper critical dimension 
is given. 
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